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It was recently realized that three-flavor effects could peculiarly modify the development of 
spectral splits induced by collective oscillations, for supernova neutrinos emitted during the cooling 
phase of a protoneutron star. We systematically explore this case, explaining how the impact of 
these three-flavor effects depends on the ordering of the neutrino masses. In inverted mass hierarchy, 
the solar mass splitting gives rise to instabilities in regions of the (anti)neutrino energy spectra that 
were otherwise stable under the leading two-flavor evolution governed by the atmospheric mass 
splitting and by the 1-3 mixing angle. As a consequence, the high-energy spectral splits found 
in the electron (anti)neutrino spectra disappear, and are transferred to other flavors. Imperfect 
adiabaticity leads to smearing of spectral swap features. In normal mass hierarchy, the three-flavor 
and the two-flavor instabilities act in the same region of the neutrino energy spectrum, leading to 
only minor departures from the two-flavor treatment. 

PACS numbers: 14.60. Pq, 97.60.Bw 



I. INTRODUCTION 

The neutrino flux from a core-collapse supernova 
(SN) is a powerful tool to probe fundamental neutrino 
properties as well as the dynamics of the explosion [l], 
1^]. The diagnostic role played by neutrinos during a 
stellar collapse is largely associated with the signatures 
imprinted on the observable SN neutrino burst by flavor 
conversions occurring deep inside the star. 

It has been understood that the paradigm of neutrino 
flavor transformation in supernovae 0], based primarily 
on the Mikheyev-Smirnov-Wolfenstein (MSW) effect 
with the ordinary matter was incomplete. New 
surprising and unexpected effects have been found to be 
important in the region close to the neutrinosphere (see 
0, i6(| for recent reviews) . Here the neutrino density is 
so high that the neutrino-neutrino interactions dominate 
the flavor evolution, producing collective oscillations. 
The most important observational consequence of v-v 
interactions is a swap of the Ve and Vf, spectra with 
the non-electron Vx and Vx spectra in certain energy 
ranges 0,lll- 

The development of the spectral swaps is strongly 
dependent on the original SN neutrino fluxes. In the 
accretion phase, i.e. t 0.5 ms after the core-bounce, 
neutrino number fluxes are expected to be ordered as 



In such a scenario, one 



flnds that for normal neutrino mass hierarchy (NH, 
Am^j^ = m| — m\ 2 > 0) collective oscillations do 
not play a significant role. For inverted hierarchy (IH, 
Am^tni < 0), the end of collective oscillations is marked 
by a complete exchange of the e and x flavors for almost 
all antineutrinos. For neutrinos, the exchange happens 
only above a characteristic energy fixed by the lepton 
number conservation, giving rise to a spectral split in 



their energy distributions pJHli 

The neutrino number fiuxes may be significantly 
different at later times, i.e. during the cooling phase. In 
Garching simulations one finds a cross-over among 
the different v spectra. As a consequence, the original 
fluxes exhibit a different ordering <I>°^ > > 
A study of this latter case was performed [l^, finding 
the occurrence of unexpected multiple spectral splits for 
both neutrinos and antineutrinos, in normal and inverted 
mass hierarchies. The rich phenomenology of the spectral 
splits, and its dependence on the original neutrino energy 
spectra was further explored in the extensive study 
performed in p7| . 

Most of the collective flavor dynamics of SN neutrinos 
can be explained in an effective two-flavor {2v) frame- 
work [1^. Collective oscihations are triggered by an 
instability in the two-flavor "iJ-sector" associated with 
the atmospheric mass-squared difference Am^^j^j and the 
mixing angle Oi^. Three-flavor {'iv) effects are due to the 
"L-sector," governed by the smaller solar mass splitting 
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TTij — m{ > and the mixing angle 612- These 



effects have been studied for neutrino fluxes typical of the 
accretion phase (l8l - l2l| . It was recently shown 22 1 that 
they are able to trigger collective flavor conversions, even 
if the mixing angle ^13 is exactly zero. However, apart 
from this initial kick, no new sizeable effect was found in 
the subsequent neutrino flavor evolution. 

Recently, the case was studied for a scenario 
relevant to the cooling phase [2^. It was found that 
in the inverted mass hierarchy, the presence of the 
solar sector can "erase" the high-energy spectral splits 
that would have occurred in the Vf, and v^. spectra 
for a 2v flavor evolution governed by the Am^^^ and 
6*13. Moreover, the final electron antineutrino energy 
spectrum exhibits a "mixed" nature, i.e. the spectral 
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swap is not complete. The phenomenological importance 
of this result is underlined by the fact that the high- 
energy spectral features arc expected to be si gnif icantly 
easier to observe at neutrino detectors. In [2j] these 
three-flavor effects are associated with an instability in 
the L sector and to the subsequent non-adiabatic flavor 
evolution driven by Am^^j. Building on this insight, we 
find it worthwhile to take a closer look at the three-flavor 
effects during the cooling phase, and to understand the 
origin and nature of the instability. 

We explain how the effects are crucially dependent 
on the neutrino mass hierarchy. For normal mass 
hierarchy, there is no fundamental difference between H 
and L sectors, since in this case both mass splittings 
are positive and the effective in-medium mixing angles 
are both small. Thus, even in a 2v set-up one would 
expect spectral splits driven by the L-sector parameters, 
and expect them to be exactly where the splits appear 
in the case of i7-sector in normal hierarchy. However 
we find that, it is not the case. In fact with typical 
SN parameters, collective oscillations fail to produce 
any significant conversions for the L-sector. This is 
so, because the L-sector has a lower natural frequency 
(wL = Aml^j2E, where E is a typical SN v energy) than 
the 77-sector and the collective interaction strength drops 
at a rate much faster than it. This does not leave enough 
time for the instability to grow and makes the collective 
evolution non-adiabatic. Therefore, spectral splits fail 
to develop. Nevertheless, the system is in an unstable 
situation. Indeed, as we will show, a small perturbation 
in the initial conditions is enough to develop the collective 
flavor conversions, producing spectral swaps also in this 
two- flavor case. 

In a realistic situation, this initial perturbation for 
the L-sector is provided by iv effects that couple this 
sector to the i/-sector. Oscillations in the i7-sector are 
communicated to the L-sector, allowing the instability to 
grow much faster. The spectral swapping still remains 
less adiabatic than in the _ff-scctor. As a result, for 
normal mass hierarchy - where both H and L-sectors 
trigger the same instability and "compete" to convert the 
high-energy v and v spectra - the iJ-sector wins. Three- 
flavor effects do not cause a significant change. The 
interesting bit happens for inverted mass hierarchy. In 
this case, the 77-sector and the L-sector have instabilities 
in different parts of the spectrum and therefore do not 
compete with each other. Instead, they "cooperate" and 
act on complementary parts of the energy spectra. The 
L-sector instability catalyzed by the i/-sector, operates 
in the high-energy region without hindrance, and causes 
an additional swap, that erases the spectral split found 
in the 2v flavor evolution. The low adiabaticity in 
the L-sector is responsible for somewhat smeared splits, 
and the effect is particularly important for splits at 
higher energies, especially for the antineutrinos. In the 
remainder of this paper, we illustrate these aspects using 
simple examples and provide a semi-analytical treatment. 

The plan of our work is as follows. In Section 2, 



we present our formalism for the SN neutrino flavor 
evolution and set up our numerical calculations, i.e. 
state our inputs for neutrino masses, mixing parameters, 
original supernova neutrino energy spectra and lumi- 
nosities at late times. In Section 3, we present the 2v 
results in the H and L sectors, in particular showing 
the lack of adiabaticity in the L-scctor and the role 
of the small perturbations in the initial conditions to 
circumvent that. In Section 4, we present a complete 
2>v calculation, showing that the i?-sector catalyzes the 
L-sector and then instabilities in both sectors develop 
- in cooperation for inverted mass hierarchy, and in 
competition for normal mass hierarchy. We provide an 
estimates of the relative adiabaticity of the low energy 
and high energy splits - explaining the mixed spectra 
observed in the antincutrino sector. Finally in Section 5, 
we comment on our results and conclude. 



II. EQUATIONS OF MOTION AND 
NUMERICAL FRAMEWORK 

Mixed neutrinos are described by matrices of density 
Pp and pp for each (anti)neutrino mode. The diagonal 
entries are the usual occupation numbers whereas the off- 
diagonal terms encode phase information. The equations 
of motion (EoMs) are [8), i24,] 

i9tPp = [Hp,pp] , (1) 

where the Hamiltonian is 

Hp ^ f7p + y+V2GF j i!S_(pq_pq)(i_Vq.Vq) , (2) 

Vp being the velocity. The matrix of the vacuum oscilla- 
tion frequencies for neutrinos is fip = (mj, m|, to|)/2|p| 
in the mass basis. For antineutrinos 17p — > — fip. The 
matter effect due to the background electron density 
is represented, in the weak-interaction basis, by = 
V2Gfnediag(l,0,0). 

In spherical symmetry the EoMs can be expressed as 
a closed set of differential equations along the radial 
direction [lH, [2^. The factor (1 — Vq • Vq) in the 
Hamiltonian, implies "multi-angle" effects for neutrinos 
moving on different trajectories 0. However, for real- 
istic supernova conditions the modifications are small, 
allowing for a single-angle approximation [l^, HB]. We 
implement this approximation by launching all neutrinos 
with 45° relative to the radial directions [l9l[25j. 

For the numerical illustrations, we take the neutrino 
mass-squared differences in vacuum to be Am^^j^^ = 
2 X 10"^ eV^ and Am^^j = 8 x 10"^ eV^ close to their 
current best-fit values |27| . The values of the mixing 
parameters relevant for SN neutrino flavor conversions 
are sin^ 6*12 0.31 and sin^ 6*13 < 0.04 The matter 
effect in the region of collective oscillations (up to a few 
100 km) can be accounted for by choosing small (matter 
suppressed) mixing angles [28j , which we take to be ^13 = 
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and considering as effective mass-square 



differences in matter Ar 



and Ato^qj 
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28], 129(L We 

ignore possible subleading CP violating effects [21[ by 
setting Sep = 0. MSW conversions typically occur after 
collective effects have ceased [lO, H^l • Their effects then 
factorize and can be included separately. Therefore, we 
neglect them in the following. 

The content of a given neutrino species i^a with energy 
E a,t a position r is given by 



Paa{E,r) 



(3) 



F{E,r) 

where F^^ is the energy distribution of Va and F is the 
sum of the energy distributions of all flavors, for neutrinos 
and antineutrinos respectively. For the three relevant SN 
V energy distributions at the neutrinosphere F^^ , we take 

Fl{E) = ^l^^,SE) , (4) 

where = Li,^/{Ei,^) is the number flux, defined in 
terms of the neutrino luminosity L^^ and the neutrino 
average energy {E^^ ) . ip^^ (E) is the normalized neutrino 
spectrum {J dE tp^^^ =1), parametrized as [ll|: 



(5) 



r(/3) {E.y^ 

where is a spectral parameter, and r(/3) is the Euler 
gamma function. The values of the parameters are model 
dependent [1, For our numerical illustrations, we 

choose 

{{E,^), (E,^), {E,J) = (12, 15, 18) MeV , (6) 

and /3 = 4, from the admissible parameter ranges [lll |. 
We take ratios of the fiuxes at late-times to be UM 



$S! : $° : $" = 0.85 : 0.75 : 1.0 



(7) 



where we have assumed equal fi and r neutrino and 
antineutrino initial fluxes. 

The strength of the neutrino-neutrino interactions can 
be parametrized as 

fio = V2GF\Fl-Fl\ , (8) 

where the fluxes are taken at the neutrinosphere radius 
i? = 10 km. Numerically, we assume /io = 7 x 10^ km~^. 
When formally deriving the single-angle approximation, 
one explicitly obtains that the radial dependence of the 
neutrino-neutrino interaction strength can be written 
as Ei 



(9) 



where the r^^ scaling comes from the geometrical flux 
dilution, and the collinearity factor 
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(10) 



arises from the (1 — cos6') structure of the neutrino- 
neutrino interaction. The asymptotic behavior for large 
r agrees with what one obtains by considering that all 
neutrinos are launched at 45° to the radial-direction [l9j . 
The known decline of the neutrino-neutrino interaction 



strength, //(r) 



for r ^ i?, is evident. 



III. MULTIPLE SPECTRAL SPLITS IN A 
TWO-FLAVOR SCENARIO 

A. 21^ ff-system 

We start our investigation with the flavor evolution 
in the iJ-sector characterized by (Am^^j^^, ^13). As 
explained in jlGj . spectral swaps can develop around 
energies E = Ec of the original neutrino spectra (except 
at £' = 0), where spectra of different flavors cross each 
other, i.e. at the crossing points 

F^SEc) - F,^ {E,) = , 

F^^ {E,) - F^JE,) = , (11) 

for neutrinos and antineutrinos, respectively. A given 
crossing point is unstable if 



d{F^^ ~ F^J/dE <0 forlH, 
d{F,^- F,J/dE >0 forNH, 



(12) 



and analogously for antineutrinos. 

In Fig. 1, we show the results of the flavor evolution 
for neutrinos (left panel) and antineutrinos (right panel) 
in inverted mass hierarchy. In the upper panels we show 
the initial and flnal and energy spectra, while in 
the lower panels we show the electron neutrino survival 
probability Pee- Following the instability conditions 
stated above, one finds that in IH the spectral swap 
would develop for neutrinos around -Ec — 13 MeV and 
for antineutrinos around Ec — 9 MeV. The development 
of a spectral swap in the middle of the energy spectra 
produces two splits in the final ly spectra. In particular, 
in the i^e final spectrum the swapped region is between 
5 MeV and 23 MeV, while for P^. is between 3 MeV and 
17 MeV. 

In Fig. 2, the corresponding results for normal hierar- 
chy are shown. In this case the only unstable crossing 
point is at Pc — > 00 in the tail of the energy spectra, 
therefore the resulting ly and D spectra exhibit only a 
high-energy swap and a single split each. In this case, 
the split is at S ~ 23 MeV for ^/g, and at ~ 17 MeV 
for De- 

From this example, we realize that in the two mass 
hierarchies the instabilities occur around different and 
well-separated crossing points. This leads to spectral 
swaps that occur in different energy ranges for the 
two mass hierarchies. In fact, these ranges are non- 
overlapping and almost complementary, i.e. the high- 
energy ends of the swaps in IH are the low-energy ends 
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FIG. 1: Flavor evolution in iJ-sector for neutrinos (left 
panels) and antineutrinos (right panels) in inverted mass 
hierarchy. Upper panels: the Ve (red) and i/^: (blue) energy- 
spectra, before (dashed curves) and after (continuous curves) 
the collective oscillations. Lower panels: survival probability 
Pee for electron (anti)neutrinos. The grey-regions represent 
the range in which the spectral swap occurs. 



FIG. 2: Flavor evolution in the 21/ H-sector for neutrinos 
(left panels) and antineutrinos (right panels) in normal mass 
hierarchy. Upper panels: the Ue (red) and i/^ (blue) energy- 
spectra, before (dashed curves) and after (continuous curves) 
the collective oscillations. Lower panels: survival probability 
Pee for electron (anti)neutrinos. The grey-regions represent 
the range in which the spectral swap occurs. 



of the swaps in NH. This implies that if we take the i/ 
spectra swapped by the conversions in IH, as an input 
for conversions in NH, these latter would swap also 
the high-energy spectrum of the electron (anti)neutrinos, 
giving the impression that the high-energy split has been 
"erased." This observation will play an important role in 
our understanding of the full three-flavor evolution. 



B. 



21/ L-system 



We now consider the i-system. Since Am'^^^ > 0, wc 
expect a behavior similar to that of the i/-system in 
normal hierarchy. However, we find that in this case no 
flavor conversion occurs. This has to be attributed to two 
reasons - insufficient growth of the instability, and lack of 
adiabaticity. The strength of the instability is given by 
the off-diagonal components in the density matrix. For a 
simple pendular system with energy E, the time-period 
for growth of off-diagonal com pon ents is a few times the 
pendular time-period Tpond 
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^pcnd 



' 2E 



(13) 



which scales logarithmically with the small in- 
mcdium mixing angle. This time period is roughly 
\J Afh'l^^/Am-^^i « 8 times larger for the L-sector, 
causing the instability to develop relatively slowly. The 
slow growth is further exacerbated by a relatively fast 
decrease in collective neutrino interaction strength fj,. 

During the spectral swapping phenomenon, the spec- 
trum near the crossing acts like an inverted pendu- 
lum [l6j . The swap sweeps through the spectrum on each 
side of the crossing, and the modes at the edge of the 
swap precess at an average oscillation frequency k [16| . 
Adiabaticity requires f23j 



dlnfi 



dr 



< K. 



(14) 



Since the collective neutrino interaction strength fi goes 
as the rate at which neutrino refractive effects are 
decreasing is \d\nij,/dr\ ~ 1/50 km^^ at r ~ 200 km, 
approximately where the swapping takes place. One finds 
K sa Arh'^^i/2E ~ 1/400 km^^ for a typical energy E ~ 
32 MeV in the region of the swap, so the adiabaticity 
condition is not met. 

Nevertheless, the L-systcm has the same instability as 
the i7-system in NH, around the crossing point Ec — 
oo of the (anti)neutrino energy spectra. Indeed, if we 
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FIG. 3: Flavor evolution in the 21/ L-sector for neutrinos 
(left panels) and antineutrinos (right panels). Upper panels; 
initial Ve (red) and Vx (blue) energy-spectra and the 
spectrum after the collective oscillations, for different values 
of an initial off-diagonal component of the density matrix 
l/Oe^. I = e X (pee + Pxx) (see the text for details) . Lower panels; 
survival probability Pee for electron (anti)neutrinos. 
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FIG. 4; Growth of the off-diagonal component |pea:| of the 
neutrino density matrix for the 2v //-system in normal 
hierarchy and for the L-system taking two different values 
of I Pea; I = e X (pee + Pxx) ■ The ueutriuo energy is taken to be 
E = 32 MeV. 



of the ofF-diagonal component is almost absent for the L- 
system with e = 10~®. Even when there are adequately 
large off-diagonal elements, as in the case of e = 10""*, 
the growth is relatively slow compared to the iJ-system, 
so that significant flavor changes can start only at large 
radius ( r ^ 300 km). Consequently, due to the violation 
of adiabaticity, collective flavor conversions do not have 
enough time to develop complete splits before the effects 
of the neutrino-neutrino interactions become negligible. 



slightly perturb the initial conditions of our v ensemble, 
e.g. put explicit off-diagonal terms in the initial neutrino 
density matrix, the instability grows easily and leads to 
significant flavor conversions.^ 

In Fig. 3, we show the effects of these artificially 
triggered flavor conversions on the electron neutrinos (left 
panels) and antineutrinos (right panels). Upper panels 
show the energy spectrum, while the lower panels display 
the neutrino survival probability. We find that it is 
sufficient to take off-diagonal seeds in the density matrix 



(15) 



with e > 10^^, to produce significant flavor conversions. 
However, relative to the iJ-sector in the normal hierarchy, 
the swap is less sharp. This is naturally expected because 
the adiabaticity does not change signiflcantly by taking 
an initial perturbation. 

In Fig. 4, we show the \pex \ component of the v density 
matrix for E = 32 MeV for three cases: (i) the i?-system 
with NH, (ii) the L-system with e = 10~^, and (iii) the L- 
system seeded with e = 10~^. We realize that the growth 



^ The speed-up of flavor instabilities under the effects of very small 
seeds in the initial conditions was already pointed out in [sil . 



IV. MULTIPLE SPECTRAL SPLITS IN A 
THREE-FLAVOR SCENARIO 

A. Speed-up of the AmgQi-instability 

Equipped with the insights of the previous section, we 
arc ready to analyze the behavior of the flavor conversions 
in the complete three-flavor scenario. We work in the 
rotated basis {ve,Vx,Vy) = {92z){ve,Vfi,VT) [l^- This 
is equivalent to take 023 = in the neutrino mixing 
matrix (which makes no difference to i^e and Df, evolution, 
if p, and r flavors are treated identically). The vacuum 
Hamiltonian is written as llSll 



n{E) 



Am 



-2 



2E 



2E 



sis C13S13 


C13S13 C?3 



(16) 



^13*12 
C12C13S12 
-C13SI2S13 



C12C13S12 

'-12 
-C12S12S13 




where Cij = cos 



6*13 = 0, the i^e 



9ij and sij = sinSy . In the limit of 
(H) and the i^e — (L) sectors are 
completely decoupled. 

When the neutrino-neutrino interaction is strong, all 
the density matrix mode pp stay pinned to each other, 
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exhibiting collective flavor conversions. Now we review 
the factorization of the H sector from other sub-leading 
contributions, as shown in For each energy mode, 

we rewrite our density matrix and Hamiltonian as 



H = H(") + hW, 



(17) 
(18) 



where the superscript (0) refers to off-components in the 
L-sector and all diagonal components, while (1) to all 
others, namely 



JO) 



and 



Pee Pex 
P*ex Pxx 

Pyy 



Pey 
pxy 

Pey Pxy ^ 



(19) 



(20) 



Analogous expressions hold for H'"^ and H'^^'. 

Note from Eq. (HH) that H(°), which contains the 
e — X block, causes oscillations due to Arn^^^, while H^^^ 
containing the e — y off-diagonal terms, gives Ato^^^- 
driven oscillations. Putting this decomposition into the 
equations of motion, one finds [l^ 



(1) 



(21) 
(22) 



Analogous equations hold for the antineutrinos. Note 
the interesting structure of the EoMs, that is an 



outcome of the commutation relations 
only on commutators [H^'^-', and 
p^^-* depends only on the cross-terms [1 
In a pure 2^ L-systcm evolution, 




depends 
while 



zero. Thus Eq. (|2ip has only the first term on r.h.s., and 
Eq. is irrelevant. Once off-diagonal components of 
p^^\ driven by develop asymmetrically in neutrinos 
and antineutrinos, the instability grows under the action 
of H^^'' (proportional to Ato^^j), i.e. the pendular time- 
scale is TVjnni 



The addition of off-diagonal 
elements to p^^^ kick-starts the process. 

In a complete system, H^^^ is non-zero and thus 
produces the off-diagonal components in p^^^ and /o'-^-' 
more quickly. Not only do initial off-diagonal terms get 
generated in the L-sector, but also those terms grow at 
a much faster rate. The growth is speeded up by the 
second term in the equation of motion Eq. (|2ip . i.e. by 



H(^) which induces oscillations Am^jj^-depcndent at the 
leading order. Therefore, it leads to a growth of the 
instability at almost the same rate as for the i/-system, 
and much faster than an isolated L-system. 

In Fig. 5, we show the diagonal components of the 
neutrino density matrix pee, Pxx, Pyy (upper panel) and 
the off-diagonal p^x , Pey (lower panel) for a given energy 
mode with _B = 32 MeV as a function of r in inverted 
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FIG. 5: Initial radial evolution of the diagonal pee, Pxx, 
Pyy (upper panel) and off-diagonal pex, Pey (lower-panel) 
components of the neutrino density matrix for inverted mass 
hierarchy, shown for a neutrino energy E — 32 MeV. 



mass hierarchy. The initial behavior is qualitatively 
similar to the one in normal hierarchy (not shown). 
The pendular oscillations of pee begin to develop at 
r ~ 35 km for the e — y sector, and proceed as pure 2v 
transitions till r ~ 50 km. Up to this point, pex has 
not evolved significantly. All off-diagonal components 
increase rapidly, but the pex starts to develop only after 
Pey , and saturates at r ~ 60 km, when e — x conversions 
start. Note that, pex grows faster than a 2i/ L-system 
(shown in Fig. 4), as predicted. The 0i3-coupling between 
the L and H sectors induces three-flavor effects in the 
neutrino conversions: The initial kick, associated with 



system. 



is necessary to trigger the instability in the L- 



B. Inverted mass hierarchy 



In Fig. 6, wc show the complete development of the 
Pee, Pxx and Pyy components of the density matrix 
in inverted mass hierarchy for neutrinos (left panels) 
and antineutrinos (right panels) for four representative 
energy modes. We observe that once the conversions 
have been started, the different energy modes for the 
three neutrino species oscillate in phase, confirming the 
collective behavior of the flavor conversions. However, 
the final fate of pee , Pxx and pyy depends on their energy. 
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FIG. 6: Three- flavor evolution in inverted mass hierarchy. 
Radial evolution of the diagonal components of the density 
matrix p for neutrinos (left panels) and antineutrinos (right 
panels) for different energy modes. 



As a result of the three-flavor effects, the v^. mixes with 
both Vx and Uy. Therefore, the t/g flavor conversions can 
be described in terms of the combined effects of the L 
and H two-neutrino systems. We find that the effects of 
the i/-sector on the v,, — Vy conversions saturate before 
the ones of the L-sector, as expected by the hierarchy 
between the two mass splittings. As we have discussed in 
Section III-A, Am^^j^^ < and Am^^j > are expected 
to process complementary parts of the neutrino energy 
spectra. Therefore, their effects do not interfere in the 
same energy range. In particular, in the part of the 
neutrino energy spectra unstable under the effects of 
Am^j^, the pee and pyy would swap, while pxx, which has 
been perturbed from its stable equilibrium configuration, 
would relax to it again. Conversely, in the part of the 
neutrino energy spectra unstable under the effects of 
Am^^j, the pee and p^x swap, while pyy comes back to its 
initial value. 

In the three-flavor space, the neutrino ensemble 
behaves like a pendulum. Once it is perturbed from its 
initial configuration, it would evolve toward its stable 
equilibrium position which may be different for different 
energy modes. In inverted mass hierarchy, the highest 
energy modes relax to the x state, while the intermediate 
ones to the y state, while lowest energy modes remain in 
the e flavor. 

For the modes represented in Fig. 6, we see that for 



FIG. 7: Three-flavor evolution in inverted mass hierarchy 
for neutrinos (left panels) and antineutrinos (right panels). 
Upper panels: Energy spectra initially (dashed curves) and 
after collective oscillations (solid curves) for (red), 
(black) and Vy (blue). Lower panels: probabilities Pee (solid 
red curve). Pet/ (dashed blue curve), P^y (dotted black curve). 
The grey bands represent the region where the spectral swaps 
occur. For antineutrinos the light grey band indicates the 
region where the spectral swap is partial (see text for details). 



E = 2.5 MeV none of the three neutrino species is 
affected by significant flavor changes. For = 7.5 MeV 
the flavor conversions produce a swap between pee and 
Pyy, while pxx comes back to its original value. At 
23 MeV, we are in the transition region between e — y 
and e — X conversions, therefore pyy and pxx are both 
partially swapped into pee- Finally, at i? = 40 MeV, pee 
and Pxx exchange their initial values, while pyy returns to 
its original value. Therefore, the behavior of the diagonal 
components of the density matrix at different energies 
would produce a single split in the v^, energy spectrum, 
since all the Ve modes at sufficiently high energy would 
swap with either Vy or Vx- Therefore, the high-energy 
spectral split at £' ~ 23 MeV observed in inverted mass 
hierarchy in the 2v evolution of Sec. III-A is washed-out 
by the three-flavor effects. In the antineutrinos (right 
panels), we find an analogous behavior. Kt E = 2.5 MeV 
Pee, Pxx and Pyy are essentially unchanged. At E = 
13 MeV the flavor conversions produce a swap between 
Pee and Pyy, while pxx comes back to its original value. 
For i? = 23 MeV, Pyy returns to its original value, 
while Pee and pxx tend to exchange their initial values. 
However, as we will discuss later, due to the the less 
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adiabatic behavior around the splitting region, the v^-Vx 
swap is not complete. Finally, for E = AQ MeV Pee and 
Pxx completely convert into each other and pyy is stable. 

In the upper panels of Fig. 7, we show the neutrino (left 
panels) and antineutrino (right panels) spectra before 
and after the complete Zv flavor conversions in inverted 
hierarchy. In the lower panel we show the corresponding 

Pee = P{Ve ^ Ve), Pex = P{Ve Vx) and Pey = P{Ve ^ 

Vy) probabilities, which can be defined approximately as 
the off-diagonal elements of the density matrix are small 
at the end of the collective evolution. Namely, these three 
probabilities read 



P{ve Va) 



raa 



P° 



rxx 



(23) 



where and / indicate the initial and final values of 
p, respectively. We first consider the neutrino flavor 
evolution. As we have already discussed, the effects 
of the L and H systems process complementary parts 
of the neutrino energy spectra. Indeed, as we can see 
from the conversion probabilities, in the energy range 



where Pp 



1, wc have Pp, 



and viceversa. In 



particular, the solar mass splitting Am^^j induces an 
electron survival probability Pee = at high-energies, 
erasing the iJ-sector high-energy split at ~ 23 MeV, as 
was pointed out in [23j . Thus, the final electron neutrino 
spectrum reads 



ioT: E<h MeV 
for > 5 MeV , 



(24) 



that occur later can make the high-energy split reappear 
in the observable electron neutrino spectra at Earth. 

Moving on to the antineutrinos, we realize that the e—x 
swap is not sharp, due to the imperfect adiabaticity in the 
L sector. This implies that the swap of De and Vx spectra 
at intermediate energies (14 MeV < < 40 MeV) is 
not complete and this explains the "mixed" Ve spectrum, 
observed also in the numerical simulations in [23]. This 
effect has to be attributed to the low adiabaticity of the 
Amg^j-induced conversions, which is particularly visible 
for antineutrinos. The adiabaticity of a spectral split 
depends on the condition in Eq. (|14p . The i>e and Vx 
spectra are closer to each other than the corresponding 
neutrino spectra. This implies that in the high-energy 
region, where the swap takes place. 



\F,^ {E) - F,^ {E)\ « \F,^ (E) - F,^ iE)\ 



(27) 



In the final phases of the swapping dynamics, the 
neutrino and antineutrino spectra evolve quite indepen- 
dently, and the precession frequencies of the two blocks 
are not governed by a common p, but by individual /i's 
proportional to the flux differences in Eq. ^7} . They 
behave essentially as two uncoupled oscillators because 
the neutrino-neutrino interaction p is now smaller than 
the frequency difference of the two blocks. One can see 
this clearly in the numerical simulations of Ref. p^ . 
shown at [11]. The frequency t~^^^ w is lower 

for antineutrinos [l6j . and thus adiabaticity tends to be 
broken more severely for antineutrinos. 



C. Normal mass hierarchy 



Conversely, the behavior of the y and x spectra is the 
same as the that of the non-electron species in 2i^ flavor 
evolution in IH and NH, respectively, i.e. 



Pf 



PeyFl + (1 



p "if" 

for E <5 MeV 
for 5 MeV < £; < 23 MeV 
for E>23 MeV , 



(25) 



and 



P 



(26) 



F^ for < 23 MeV 
F,f for E > 23 MeV . 



Indeed; the Vx and Vy spectra will be affected respectively 
only by the L or the H sector, since x ^ y conversions 
are strongly suppressed as evident from the vacuum 
Hamiltonian in Eq. One should note that while 

the high-energy spectral split is no longer present in 
the Ve final energy spectrum, Vy and Vx final spectra 
still present this high-energy feature. Therefore, MSW 
effects in supernova, vacuum mixing and Earth effects 



In Fig. 8, wc show the complete development of the pee, 
Pxx and Pyy components of the density matrix in normal 
hierarchy for neutrinos (left panels) and for antineutrinos 
(right panels) for different energy modes. For both 
neutrinos and antineutrinos, the swapping dynamics 
occurs mainly between the pee and the pyy^ while pxx 
experiences only nutations around its initial equilibrium 
value, before relaxing completely to it when the flavor 
conversions are saturated. Therefore, the flavor evolution 
is close to the 2v H-case in normal hierarchy discussed in 
Sec. III-A. Indeed, both Am^^jj^ and Am^^j are positive. 
Therefore, the L-sector would behave as a replica of the 
ff-sector, but with a smaller mass splitting. In this 
condition, both H and L sectors process the same regions 
of the electron neutrino energy spectra. However, the 
hierarchy between the two mass splittings produces the 
dominance of the Ue — Vy swaps, while conversion effects 
in the Ve — Vx sector remain inhibited by adiabaticity 
violation. The only region where the L-instability can 
compete with Fl is very close to the split. In particular, 
for low-energy neutrino modes {E = 2.5, 18 MeV in 
Fig. 8) the pee comes back to its initial value, while 
at higher energies (E = 23, 40 MeV) exchanges its 
initial value with pyy. For the antineutrinos, Pee remains 
at its initial value for the modes at E = 2.5 MeV. 
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FIG. 8: Three- flavor evolution in normal mass hierarchy. 
Radial evolution of the diagonal components of the density 
matrix p for neutrinos (left panels) and antineutrinos (right 
panels) for different energy modes. 



FIG. 9: Three-flavor evolution in normal mass hierarchy 
for neutrinos (left panels) and antineutrinos (right panels). 
Upper panels: Energy spectra initially (dashed curves) and 
after collective oscillations (solid) for Ue (red), Ux (black) 
and Vy (blue). Lower panels: corresponding probabilities 
Pee (solid red curve), Pey (dashed blue curve), Pey (dotted 
black curve). The grey bands represent the region where the 
spectral swaps occur (see text for details). 



At E = 18 MeV, we are around the splitting region 
where the adiabaticity in the H-systeui is more severely 
violated, and the e-y conversions are not complete. 
Under these conditions, the L instability can also play 
a role, producing a weak swap in p^x- Finally, at 
_E = 40 MeV, conversions again occur only between e 
and y states, producing a complete swap of the initial 
and pyy values. 

In Fig. 9, we show the neutrinos (left panels) and 
antineutrino (right panels) spectra before and after the 
complete 3i> flavor conversions in NH. In the lower 
panels we represent the corresponding Pee, Pex and P^y 
probabilities. Once more, we realize that the flavor 
evolution can be mostly described in terms of two- 
flavor i/g — Vy transitions, while the role of v,, — 
conversions is subleading. The spectral splitting features 
are remarkably similar to the ones found for the 2v H- 
system in normal hierarchy (Sec. Ill- A). In particular, the 
i^e spectrum presents a single split at i? ~ 23 MeV and i>e 
spectrum around i? ~ 17 McV. Only close to the splitting 
regions, we find subleading effects associated with i/g — Vx 
conversions. 



V. CONCLUSIONS 

Collective neutrino flavor conversions in supernovae, 
associated with neutrino-neutrino interactions, have been 
recognized to induce peculiar spectral swaps among the 
different neutrino species. The development of these 
features is associated with instabilities in the flavor 
space. In particular, these instabilities would develop 
around the crossing points of the original SN neutrino 
spectra. Then, the neutrino mass hierarchy determines 
if a crossing point is unstable under the effects of the 
collective oscillations. A particularly intriguing case is 
the one in which the original SN neutrino fluxes exhibit 
an ordering with $° $2 ^ ^1 i possible during the 
coo ling phase. In this case, the two-flavor study realized 
by [l6|, found the occurrence of multiple spectral splits 
for both neutrinos and antineutrinos, depending on the 
neutrino mass hierarchy. A recent numerical exploration 
of this case performed in [2^ , has found that when three- 
flavor effects are taken into account in inverted mass 
hierarchy, the high-energy spectral swaps observed in the 
21/ evolution are erased by effects related to Amg^j. 

Motivated by this intriguing result, in our paper we 
have performed a detailed study of the three-flavor 
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effects in the collective oscillations for supernova neutrino 
spectra typical of the cooling phase. We have found 
that the effects of ATOg^j in the three-flavor evolution 
are important only in inverted mass hierarchy. In this 
case, the presence of Am^^j gives rise to instabilities in 
regions of the neutrino energy spectra that were stable 
under the two-flavor evolution governed by Am^^j^^ and 
6*13. Therefore, the combinations of these two different 
instabilities would produce a wash-out of the high-energy 
splitting features in the and v^, spectra. Conversely, in 
normal mass hierarchy the three-flavor instabilities and 
the two-flavor one act in the same regions of the neutrino 
energy spectrum, leading only to minor departures from 
the two-flavor evolution. Essentially, the system behaves 
like a pendulum in Zv flavor space. It can topple 
towards either the Vy state or the Vx state. In inverted 
hierarchy, when the H and L instabilities are in different 
regions of energy, the pendulum topples towards Vy for 
the i7-instability, and towards for the L-instability. 
In normal hierarchy, when the instabilities are in the 
same region of energy, the pendulum topples towards 
Vy^ as the L-instability is relatively non-adiabatic. As 
a consequence, in inverted mass hierarchy the electron 
(anti)neutrino spectrum at the end of the collective 
oscillations would present only a very low-energy {E < 
5 MeV) splitting feature, being completely swapped to 
the original non-electron spectra at higher energies. 

We wish to emphasize that the high-energy split- 
ting features may survive in the observable electron 
(anti)ncutrino spectrum at Earth, even in inverted 
hierarchy. Indeed, MSW matter effects in SN, vacuum 
mixing and Earth effects would further mix the and 
Vx spectra. The non-electron spectrum at the end of the 
collective oscillations, still contains high-energy splitting 
features, since for the non-electron species the collective 
flavor conversions have occurred essentially as in the two- 
flavor case. Therefore, especially for neutrinos, which 
have sharper spectral swaps, the electron neutrino signal 



at Earth could still present observable splitting features 
at high energy. Also, due to the lower neutrino luminosity 
at sufficiently late times, the initial collective interaction 
strength //q can be somewhat lower than is assumed 
in Ref.[23[ and this work. We found that in certain 
regions of the spectral parameter space, reducing /ig by 
a factor of 10 makes the three-flavor effect disappear due 
to a stronger adiabaticity violation in the L sector. In 
principle, this effect could produce interesting signatures 
in the time evolution of the SN neutrino signal. 

In conclusion, the non-linear equations that govern 
the flavor evolution of neutrinos emitted during a stellar 
collapse are a continuous source of surprises and new 
effects. During this last year, dramatic changes have 
occurred in the picture consolidated after the initial 
exploration of collective supernova neutrino oscillations. 
The discovery of this new three-flavor effect is the most 
recent of these changes. After our study, it appears that 
its impact on the collective neutrino flavor conversions is 
conceptually and quantitatively well under control. 
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